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Type la Supernovae (SNe la) form an observationally uniform class of stellar explosions, in 
that more luminous objects have smaller decline-rates 1 . This one-parameter behavior al- 
lows SNe la to be calibrated as cosmological 'standard candles', and led to the discovery 
of an accelerating Universe 2 3 . Recent investigations, however, have revealed that the true 
nature of SNe la is more complicated. Theoretically, it has been suggested 4 8 that the ini- 
tial thermonuclear sparks are ignited at an offset from the centre of the white-dwarf (WD) 
progenitor, possibly as a result of convection before the explosion 4 . Observationally, the di- 
versity seen in the spectral evolution of SNe la beyond the luminosity decline-rate relation is 
an unresolved issue 9 ' 10 . Here we report that the spectral diversity is a consequence of ran- 
dom directions from which an asymmetric explosion is viewed. Our findings suggest that the 
spectral evolution diversity is no longer a concern in using SNe la as cosmological standard 
candles. Furthermore, this indicates that ignition at an offset from the centre of is a generic 
feature of SNe la. 
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When a carbon-oxygen WD reaches a critical limit known as the Chandrasekhar mass (~ 
1.38 M ), its central density and temperature increase to a point where a thermonuclear runaway 
is initiated. The thermonuclear sparks give birth to a subsonic deflagration flame, which at some 
point may make a transition to a supersonic detonation wave that leads to the complete disruption 
of the WD 1112 . The thermalization of 7-rays produced from the decay of freshly synthesized 
radioactive 56 Ni powers the transient source, known as a SN la 13 ' 14 . The relationship between the 
luminosity and the decline-rate parameter (Am 15 (B), which is the difference between the 5-band 
magnitude at peak and that measured 15 days later) is interpreted to be linked to the amount of 
newly synthesized 56 Ni (refs. 15, 16). 

SNe la displaying a nearly identical photometric evolution can exhibit appreciably different 
expansion velocity gradients (v S i) as inferred from the Si II A6355 absorption feature 1017 . More 
specifically, objects that show v S i > 70 km s _1 day -1 are placed into the high-velocity gradient 
(HVG) group, while those that show smaller gradients are placed in the low-velocity gradient 
(LVG) group. For normal SNe la 18 , which are the predominant population of the total SN la 
sample and the main focus of this Letter, v Si is not correlated with Am 15 (5) (ref. 10; Fig. la, lb), 
thus raising a nagging concern regarding the 'one-parameter' description. 

Late phase nebular spectra can be used to trace the distribution of the inner ejecta 19 . Begin- 
ning roughly half a year after explosion, as the ejecta expands, its density decreases to the point 
where photons freely escape. Photons originating from the near/far side of the ejecta are detected 
at a shorter/longer (blue-shifted/red-shifted) wavelength because of Doppler shifts. For SNe la, 
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emission lines related to [Fe II] A7155 and [Ni II] A7378 are particularly useful, as they are formed 
in the ashes of the deflagration flame 19 . These lines show diversity in their central wavelengths 
- blue-shifted in some SNe la and red-shifted in others (see Fig. lc) - which provides evidence 
that the deflagration ashes, therefore the initial sparks, are on average located off-centre. The 
wavelength shift can be translated to a line-of-sight velocity (v nch ) of the deflagration ashes. 

Figure 2 shows a comparison between v S i and v nc b for 20 SNe la. Details regarding the data 
are provided in SI §1. Although the diversities in these observables were discovered independently, 
Fig. 2 clearly shows that they are connected. Omitting the peculiar SN 2004dt (Fig. 2 Legend; SI 
§1), all 6 HVG SNe show v nch > km s _1 (i.e., red-shifts), which means that these events are 
viewed from the direction opposite to the off-centre initial sparks. The 11 LVGs display a wider 
distribution in v ne b space, but are concentrated to negative values (i.e., blue-shifted), indicating that 
these events are preferentially viewed from the direction of the initial sparks. If HVG and LVG 
SNe were intrinsically distributed homogeneously as a function of v nc b, the probability that just by 
chance (as a statistical fluctuation) all HVG SNe show v nc b > km s _1 is merely 0.4%; it is thus 
quite unlikely. Indeed, the probability that by chance 6 HVG SNe are among the 7 SNe showing 
the largest red-shift in t> ncb in our sample of 17 SNe is only 0.06%. 

This finding strongly indicates that HVG and LVG SNe do not have intrinsic differences, but 
that this diversity arises solely from a viewing angle effect. Figure 3 shows a schematic picture. 
If viewed from the direction of the off-centre initial sparks, an SN la appears as an LVG event at 
early phases and shows blue-shifts in the late-time emission-lines. If viewed from the opposite 
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direction, it appears as an HVG event, and shows red-shifts at late phases. 

The number of HVG SNe is - 35% 10 of the total number of HVG and LVG SNe. To explain 
this, the angle to the observer at which an SN changes its appearance from an LVG to an HVG 
is ~ 105 — 110°, measured relative to the direction between the centre and the initial sparks. 
The velocity shift of 3, 500 km s _1 in the distribution of the deflagration ashes, as derived for the 
normal SN la 2003hv through a detailed spectral analysis 19 , corresponds to t> ncb ~ 1000 km s _1 if 
viewed from this transition angle. The configuration then predicts that all LVG SNe should show 
-3, 500 < v nch < 1000 km s~\ while all HVG SNe should be in the range 1000 < v nch < 3, 500 
km s" 1 . These ranges are shown as arrows in Fig. 2, and provide a good match to the observations. 

Figure 4a shows an example of a hydrodynamic model in which the thermonuclear sparks 
were ignited off-centre in a Chandrasekhar-mass WD 6 (an alternative way of introducing global 
asymmetries is double detonations in sub-Chandrasekhar-mass WDs 20 ). Although this model has 
not been fine-tuned to reproduce the present finding, it does have the required generic features. The 
density distribution is shallow and extends to high velocity in the direction opposite to the initial 
sparks (Fig. 4b). Initially, the photosphere is at high velocity if viewed from this direction, as the 
region at the outer, highest velocities is still opaque. Later on, the photosphere recedes inwards 
faster in this opposite direction, owing to the shallower density gradient. As a result, the SN looks 
like an LVG if viewed from the offset direction, but like an HVG SN from the opposite direction 
(Fig. 4c), as in our proposed picture (Fig. 3). 

Our finding provides not only strong support for the asymmetric explosion as a generic fea- 
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ture, but also constraints on the still-debated deflagration-to-detonation transition. In this particular 
simulation, the change in appearance (as an HVG or an LVG SN) takes place rather abruptly around 
the viewing direction of ~ 140°. Owing to the offset ignition, the deflagration flame propagates 
asymmetrically and forms an off-centre, shell-like region of high density deflagration ash. The det- 
onation is ignited at an offset following the deflagration, but tries to expand almost isotropically. 
However, the angle between 0° and 140° is covered by the deflagration ash, into which the strong 
detonation wave (fueled by the unburned material near the centre of the WD) cannot penetrate. On 
the other hand, in the 140 — 180° direction, the detonation can expand freely, creating a shallow 
density distribution. The 'abrupt' change in appearance, as inferred by the observational data, is 
therefore a direct consequence of the offset models, controlled by the distribution of the deflagra- 
tion ash. The 'opening angle' of the transition is on the other hand dependent on the details of 
the explosion. To accurately model this according to our finding (~ 105 — 110° for the typical 
transition angle), either a smaller offset of the initial deflagration sparks or an earlier deflagration- 
to-detonation transition would be necessary. Such changes are also required to produce the typical 
velocity shift of ~ 3, 500 km s _1 in the distribution of the deflagration ash. 

Our proposed model unifies into a single scheme recent advances in both theoretical and 
observational studies of SNe la - and it does not conflict with other results produced by spectral 
tomography 16,21 or polarization measurements 22 (SI §2). Our interpretation suggests that two SNe 
with very similar light curve evolution may not necessarily produce exactly the same amount of 
56 Ni. They might in fact show somewhat different light curve evolution if viewed from the same 
direction (relative to the offset between the centre and the initial sparks), but the light curves look 
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the same to an observer as Am 15 (5) can change with viewing direction 5 (SI §3). This situation is 
a natural consequence of observing a number of SNe la with a large variation in 56 Ni production 
and in the viewing angles. Our finding regarding the explosion mechanism will lead to quantitative 
evaluation on the contribution of this random effect to the observed scatter in the SN la luminosities 
beyond the one-parameter description 5 , as compared to other systematic effects, such as the stellar 
environment 23 . 
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Figure 1 Comparison between HVG SN la 2002bo and LVG SN la 1998bu. The 
decline-rate parameter Am 15 (B) is 1.16 and 1.06 mag for SNe 2002bo and 1998bu, re- 
spectively, a. The S-band light curves 17,24 . The magnitudes for SN 1998bu have been 
artificially shifted in the y-direction for presentation, b. Si II A6355 at different epochs 17,24 
(in days with respect to B-band maximum). SN 2002bo had initially a larger absorption 
velocity than SN 1 998bu, but later its velocity approached that of SN 1 998bu. The velocity 
evolved quicker and the velocity gradient (v Si ) is larger for SN 2002bo than for SN 1 998bu. 
c. [Fe II] A7155 and [Ni II] A7378 in late-time spectra 21,25 . The horizontal axis denotes the 
line velocity measured from the rest wavelength of [Ni II] A7378. The rest wavelengths of 
[Fe II] A7155 and [Ni II] A7378 are marked by dashed lines. These lines are red-shifted 
in SN 2002bo while blue-shifted in SN 1998bu. The wavelength shift indicates the line- 
of-sight velocity of the deflagration ashes {v neh ) {v neh < km s \ i.e., the blue-shift, if the 
material is moving toward us). These are the strongest lines among those emitted from 
the deflagration ash according to the previous analysis of late-time emission lines 19 . In- 
deed, there are stronger lines which do not show Doppler shifts, e.g., [Fe III] A4701 ; they 
however do not trace the distribution of the deflagration ash 19 (see also SI §1) and are 
thus not useful in the present study. 
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Figure 2 Relations between the features in early- and late-phases, a. Early-phase 
velocity gradient (v Si : vertical axis) as compared to late-phase emission-line shift velocity 
(f neb : horizontal axis) for 20 SNe la. The errors are for la in fitting the velocity evolution 
for v Si , while for v nch the errors are from differences in measurement between different 
emission lines (see Sup §1). LVG SNe and HVG SNe are shown by green squares and 
magenta circles, respectively SNe 1 986G and 2007on are classified as 'faint and fast- 
declining' SNe la 26-28 (black open circles). SN 2004dt (orange square) is an HVG SN 
according to the value of i; Si , but displayed peculiarities in the late-time spectrum 29 (SI §1) 
and in polarization measurements 22 (SI §2), and the value for t; Si is exceptionally large 
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as compared to other HVG SNe. These suggest that SN 2004dt is an outlier and the 
origin of its large v Si is probably different from that of other SNe la. The two arrows 
on top indicate the regions where HVG and LVG SNe are expected, based on a simple 
kinematic interpretation (see main text), b. Number distribution of 20 SNe as a function 
of v S i- White and orange areas are for faint SNe and SN 2004dt. The remaining SNe 
are marked depending on whether they show a blue-shift {v neh < km s _1 : blue area) or 
red-shift {v nch > km s _1 : red area) in their late-time spectra, c. Number distribution of 
20 SNe as a function of v neh . 
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Figure 3 A schematic picture of the structure of SN la ejecta. This configuration 
simultaneously explains the relative fractions of HVG and LVG SNe la, and the relation 
between v Si and v nch . It is also consistent with the diversity in v nch (ref. 19). The ashes 
of the initial deflagration sparks are shifted with respect to the center of the SN ejecta by 
~ 3,500 km s 1 (yellow: Although expressed by a spherical region for presentation, it may 
well have an amorphous shape owing to the hydrodynamic instability in the deflagration 
flame 7 ). This region is rich in stable 58 Ni with a small amount of radioactive 56 Ni, and emits 
[Fe II] A7155 and [Ni II] A7378 at late phases 19 (SI §1). The outer region is the later deto- 
nation ash responsible for the early-phase Si II A6355 absorption. This region is roughly 
spherically distributed (blue), but extends to the outer region in the direction opposite to 



16 



the deflagration ashes (cyan). Although the detonation can produce 56 Ni which decays 
into 56 Fe, it has been argued that these regions (blue and cyan) are not main contribu- 
tors to [Fe II] A7155 and [Ni II] A7378 at late phases 19 (SI §1). A putative observer would 
view this explosion as an LVG or HVG SN, depending on the direction of the observer, as 
divided by a specific angle (red), which is ~ 105 - 110°. This angle is derived from the 
relative numbers of HVG and LVG SNe (see main text); the fraction of HVG SNe to the 
total number of HVG and LVG SNe ~ 35% (ref. 10), as was also supported by a larger 
sample 30 with more than 100 SNe (using the observed trend that the HVG SNe show 
higher velocities than LVG objects at early times). 
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Figure 4 Expectations from a hydrodynamic explosion model, a. Cross section of the 
density distribution of an offset ignition model 6 (similar to a model in ref. 5). It is shown 
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at 10 seconds after the ignition, when the homologous expansion is already reached. In 
this model, the deflagration sparks were ignited offset, within an opening angle of 45° 
and between - 180 km from the centre of a WD whose radius is ~ 2,000 km. The 
deflagration products are distributed in the high-density offset shell, which covers - 140° 
in this particular model. This resulting velocity shift is ~ 8,000 km s 1 in this model, 
which is larger than the observational requirement (~ 3, 500 km s 1 ). b. Radial density 
distribution of the same model for several directions (where the angle is measured from 
the direction of the offset sparks). The magenta lines show the density distribution for the 
angle in 100 - 180°, roughly corresponding to the putative directions for which a SN looks 
like an HVG SN inferred from the observational data (Fig. 3). Also shown is the density 
distribution multiplied by the mass fraction of Intermediate mass elements, which gives a 
rough distribution of Si, for two directions (20 - 40° and 140 - 160°). It is seen that the 
distribution of Si roughly follows the density distribution, c. Model photospheric velocity 
evolution for several directions (lines) as compared to the observed Si II A6355 absorption 
line evolution 10 (filled and open symbols for HVG and LVG SNe, respectively; the phase 
in days with respect to fi-band maximum). The position of the photosphere is estimated 
by integrating the optical depth with a constant opacity along each direction. Using the Si 
distribution for the velocity estimate provides a similar result. 
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Supplementary Information 



1 Supernova Sample and Notes for Individual SNe 

Supplementary Table 1 summarizes the data of the SNe used in the present study The 
values of the velocity gradients (v S i) are obtained from the literature when available. They 
are mostly drawn from a previous compilation 10 (see Supplementary Tab. 1 for the sources 
of the observations), with values for additional SNe obtained from the references listed 
in Supplementary Table 1. For SN 1998aq, we have measured v Si from the published 
spectra 38 . 

The line-of-sight velocity of the deflagration ash is derived as follows. Late-time 
spectra of SNe la show various forbidden lines from Fe-peak elements. These lines can 
be divided into two groups. The first group requires intense heating from radioactive 56 Ni 
56 Co -> 56 Fe decays and low material density (e.g., [Fe III] A4701), while the second 
group requires a low heating rate and high material density (e.g., [Fe II] A7155 and [Ni 
II] A7378). Considering general properties of the emission process 19 , the former lines are 
argued to be preferentially emitted from the detonation ash (because it is at relatively low 
density with a large amount of 56 Ni), while the latter lines are formed in the deflagration ash 
(because it is at high density with a small amount of 56 Ni). These two groups of lines show 
mutually different properties in observations, strengthening the interpretation that they are 
emitted from different regions: The 'deflagration ash'-lines show the diversity in Doppler 
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shift (not only in [Fe II] A7155 and [Ni II] A7378 61 ' 62 ), while the 'detonation ash'-lines show 
virtually no Doppler shifts. This property led to the conclusion that the deflagration ash is 
on average located offset, while the detonation ash is distributed roughly spherically 19 . 

For v ne b (the line-of-sight velocity of the deflagration ash), we have therefore mea- 
sured the wavelength shift in [Ni II] A7378 and in [Fe II] A7155, and defined v neh as the 
mean value of them (see Supplementary Tab. 1 for the source of the observational data). 
The error bars are taken to be the difference in the measurements for the two lines. In 
all SNe except for SN 2004dt, the difference in these two measurements is at most 600 
km s~ x . For SNe la in which either of [Fe II] A7155 or [Ni II] A7378 is weak and uniden- 
tified, we measure v nch only from the other single line, with a conservative error of 600 
km s _1 . For SNe 2007on, 2007sr, and 2009ab, we have measured i) Si and v nch from our 
own spectra taken at the Gemini South, the Magellan, and the ESO (La Silla, Paranal) 
telescopes. The late-time spectrum of SNe 1997bp has also been obtained by us and 
used to measure v neh . The 20 SNe in Supplementary Table 1 are all the objects we have 
found for which both v Si and v neb are reliably available to date. Also shown in the table are 
the light-curve decline-rate parameters Ami 5 (5). 

Supplementary Table 1 also lists the epoch (after S-band maximum) of the late-time 
spectra from which v nch is measured. For most of the SNe la, it is at least 200 days 
after fi-band maximum, and thus errors caused by the blending of additional lines to this 
feature (e.g., permitted lines emitted at relatively early phases) in measuring v nch can be 
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avoided. The features at ~ 7, 000 - 7, 500 A, which we interpret to be dominated by [Fe 
II] A7155 and [Ni II] A7378, do not show significant evolution at > 100 - 150 days after 
5-band maximum 19 . 

It has been shown that normal SNe la 18 (which constitute > 80% of the whole SN la 
population), faint 1991-bg like SNe la 26 , and bright 1991T-like SNe la 6364 show differ- 
ent properties in i; S i 10 . Normal SNe la show a diversity in v su which is apparently not 
correlated with Am 15 (B). Faint 1 991 bg-like SNe la, which are characterized by a rapid 
fading (i.e., large Ami 5 (B)), always show large v Si . Bright 1 991 T-like SNe la with small 
Am 15 (B) always show small v Si . In deriving t> ncb , we have noticed that the identification of 
[Fe II] A7155 and [Ni II] A7378 is robust for normal SNe la, but tends to be a subject of 
possible misidentification for the other subclasses (then such SNe la are omitted from our 
analysis): Faint 1991 bg-like SNe la tend to show probable [Ca II] A7291, 7324 27 . Bright 
1991 T-like SNe la tend to show a broad single peak, possibly indicating that [Fe II] A7155 
and [Ni II] A7378 are broad and mutually blended. These indicate that the ejecta structure 
(density, temperature, and composition) of faint and bright SNe la is intrinsically different 
from normal ones. This is a reason why our sample is mainly composed of normal SNe la. 
Luckily, omitting a large fraction of faint/bright SNe la is not important for our present anal- 
ysis, since we are interested in the spectral diversity beyond the one-parameter Am 15 (B) 
description, which is a problem only in normal SNe la. 

For normal SNe la, we categorize HVG SNe and LVG SNe according to v S u with the 
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division line at 70 km s 1 day -1 . According to its v Si , SN 2004dt is an HVG. However, 
its peculiar observational features suggest that it is an outlier, and the origins of its high 
velocity gradient and the negative v nch are likely different than the one for other HVGs 
(see Fig. 2 caption). One of the peculiar features of SN 2004dt appears in its late-time 
spectra. Supplementary Fig. 1 shows a late-time spectrum of SN 2004dt compared to 
the prototypical faint SN 1 991 bg and to the normal HVG SN 2007sr. The spectrum of 
SN 2004dt provides a good match to that of the faint SN 1 991 bg. Both the intensity and 
the width of the lines at ~ 7, 000 - 7, 500 A are similar for these two SNe. In normal 
HVG SN 2007sr, the features at 7, 000 - 7, 500 A are much fainter than for SNe 2004dt and 
1 991 bg. This indicates that in the case of SN 2004dt the feature is probably contaminated 
by other lines, most likely [Ca II] AA7291 , 7324 (note that SN 2004dt shows a larger error 
in v nch than the other SNe la). In addition, the strongest lines at late phases, i.e., the 
[Fe III] blend at ~ 4, 700 A and the [Fe II] and [Fe III] blend at ~ 5, 250 A have similar ratios 
in SN 2004dt and 1 991 bg; these ratios have been noticed to be peculiar 29 . The normal 
HVG SN 2007sr has broader lines, and different line ratios. These similarities between 
the 'normal' SN 2004dt and the faint SN 1 991 bg at late phases, which have not been 
noticed previously, suggest that these explosions might be closely related to one another. 
SN 2004dt may represent a new class of SNe la 65 which have SN 1991 bg-like features in 
the late-time spectrum, but are more energetic and brighter. 
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Supplementary Table 1: Supernovae Sample 



SN 


Am 15 {B) 


i>Si 

(km s _1 day -1 ) 


Vncb 

(km s- 1 ) 


Epoch 
(day) 


Class 8 ' 24 


References 


1986G 


1.81 ±0.07 


64 ±5 


-680 ±50 


257 


HVG/Faint 


1,31,32 


1990N 


1.07 ±0.05 


41 ±5 


-1130 ±600 


280 


LVG 


1, 10,33,34 


1994D 


1.32 ±0.05 


39 ±5 


-1440 ±390 


306 


LVG 


1, 34, 35 


1997bp 


0.97 ±0.2 


106 ±7 


2650 ± 150 


300 


HVG 


10, 36 


1998aq 


1.12 ±0.05 


35 ±4 


-1100 ±50 


241 


LVG 


37,38 


1998bu 


1.06 ±0.05 


10 ±5 


-990 ± 150 


329 


LVG 


1, 10, 24, 25,39 


2000cx 


0.93 ±0.04 


2±4 


-320 ±600 


147 


LVG/peculiar 


40-43 


200 lei 


1.13 ±0.04 


31 ±5 


1070 ± 230 


398 


LVG 


44-46 


2002bo 


1.16 ±0.06 


110 ± 7 


1900 ± 220 


368 


HVG 


17,21 


2002dj 


1.08 ±0.05 


86 ±6 


2100 ±420 


275 


HVG 


47 


2002er 


1.32 ±0.03 


92 ±5 


890 ± 600 


216 


HVG 


48,49 


2003du 


1.02 ±0.05 


31 ±5 


-370 ±50 


377 


LVG 


50,51 


2003hv 


1.61 ±0.02 


41 ±6 


-2750 ± 440 


320 


LVG 


52 


2004dt 


1.21 ±0.05 


160 ± 15 


-1910 ±1070 


152 


HVG/peculiar 


29 


2004eo 


1.45 ±0.04 


45 ±4 


0±600 


228 


LVG 


53,54 


2005cf 


1.12 ±0.03 


35 ±5 


420 ± 600 


267 


LVG 


55-58 


2006X 


1.31 ±0.05 


123 ± 10 


1280 ± 20 


277 


HVG 


59 


2007on 


1.62 ±0.01 


85 ±6 


0±600 


356 


HVG/Faint 


This work 


2007sr 


0.93 ±0.02 


80 ± 15 


1910 ±400 


256 


HVG 


60 


2009ab 


1.20 ±0.02 


36 ± 13 


-410 ±410 


278 


LVG 


This work 
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2 Other Observational Constraints 

In abundance 'tomography' 16 ' 21,66 " 68 , a temporal sequence of spectra of individual SNe la 
are used to infer the distribution of different elements through the SN ejecta, assuming the 
density structure of a spherically symmetric explosion model 13 . From this type of analysis, 
it has been indicated that the abundance distribution is generally a function of Am 15 (B). 
The difference between the HVG and LVG SNe, not related to Am 15 (B), is mainly on the 
extent of the Si-rich layer 16 and on the photospheric velocity 66 , which are explained by our 
proposed scenario (main text). 

On the other hand, the spatial extent of the 56 Ni-rich region does not seem to be 
dependent on whether it is a HVG or a LVG SN 16 . This could provide a constraint on 
the ejecta asymmetry. In the offset explosion model, the spatial extent of the 56 Ni region, 
as well as the density structure at the outer edge of that region, are not sensitive to the 
direction, despite the initially large asymmetry in the ignition 6 . This stems from the nature 
of the propagation of the detonation wave as described in the main text; unlike the defla- 
gration flame, the detonation tries to expand isotropically, producing roughly spherically 
distributed 56 Ni. This region is not sensitively affected by the existence of the deflagration 
ashes, which is essential in determining the structure of Si-rich region. As a result, the 
spatial extent of the 56 Ni-rich region is mainly controlled by the different amounts of 56 Ni 
produced in the explosion, and the viewing angle dependence could add some diversity 
at most as a secondary effect 6 ; this is consistent with the observational indications 16 . 
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The asymmetric distribution of the outermost layer may imprint its signature in po- 
larization measurements, which may be linked to the velocity gradient 69 . The polariza- 
tion of the Si II line is correlated with Ami 5 (5) 22 , but only for LVG SNe (Supplementary 
Fig. 2). HVG SNe generally show larger polarization than LVG SNe 7071 but they clearly 
do not follow this trend. A global one-sided asymmetry as in the present interpretation 
would produce relatively low continuum polarization and relatively high line polarization at 
Si II A6355 72 . In our proposed scenario the global asymmetry is of smaller degree than in 
an extremely asymmetric model producing Si II polarization of ~ 1% 72 , thus this is likely 
not a major contributor to the observed polarization. Alternatively it has been suggested, 
based on the correlation between the Si II polarization and Am l5 (B), that the observed 
Si II polarization could be a measure of the thickness of the outer layer above the 56 Ni-rich 
region, in which the local inhomogeneity, e.g., a few relatively dense blobs, is assumed to 
be a source of polarization. In this interpretation, the large Si II polarization in HVG SNe 
could be a consequence of an extended outer layer in the direction opposite to the initial 
sparks. 
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3 Viewing-Angle Effect on the Light Curve 

It has been suggested that if the ejecta are asymmetric, Am 15 (B) is dependent on the 
direction to the observer 5 . Supplementary Fig. 3 shows the comparison between Am 15 (B) 
and v nch for SNe la. 

In the same figure, the possible effect of the viewing angle on the light curve shape, 
and the implication on the scatter in the luminosity calibration beyond the one-parameter 
description, are schematically illustrated. For example, the LVG SN 1998bu and the HVG 
SN 2002bo had similar Am 15 (5) values, and this could be interpreted as follows. 

These two SNe would indeed have intrinsically different amounts of M( 56 Ni) and 
accordingly different luminosity. Because of different M( 56 Ni), Am 15 (B) would be different 
if viewed from the same direction. Assuming that (1) Am 15 (B) would be smaller and 
the luminosity is thus larger for LVG SN 1998bu than for HVG SN 2002bo, for a putative 
observer at the same direction, and that (2) Am 15 (B) appears larger for an observer 
closer to the offset direction, as a viewing angle effect. Then, these two SNe would 
show similar Am 15 (B) as observed, despite the intrinsic difference in luminosity, since the 
viewing direction is closer to the offset direction for SN 1998bu. 

This effect would produce a scatter in the luminosity of SNe la, around the standard 
luminosity determined by M( 56 Ni). Pairs of SNe with similar Am 15 (B), as a result of 
the viewing angle originating from SNe with different amounts of M( 56 Ni), are naturally 
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expected if we observe a large number of SNe la (Sup. Fig. 3) for the following two 
reasons: (1) SNe la do indeed show large variations in M( 56 Ni), and (2) the viewing 
angles should display large variations. 

A question for this interpretation is whether any indication of such an effect is seen 
in the data. Supplementary Fig. 3 shows that there is no clear (but perhaps a marginal) 
correlation between Ami 5 (B) and v nch . According to the prediction 5 of the viewing-angle 
effect on Am 15 (B) for models similar to the one shown in the main text, the observed 
Am 15 (B) could vary by ~ 0.4 mag depending on the direction to the observer (schemat- 
ically shown in Sup. Fig. 3). This is smaller than the intrinsic variation in Am lb (B) for 
different M( 56 Ni), and thus such an effect is difficult to notice in Supplementary Fig. 3, 
as is consistent with the low correlation in the present data. Any marginal correlation be- 
tween Amis (B) and v nch may already hint that such an effect is indeed there, but a larger 
number of SNe la is necessary to test this possibility with statistical significance. 
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4 Supplementary Figure 1 
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Spectrum of SN 2004dt 29 at 120 days after B-band maximum, compared to those of 
the faint SN 1991 bg 27 and of the normal HVG SN 2007sr at similar epochs. 
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5 Supplementary Figure 2 
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Si II A6355 line polarization for SNe la, as a function of the decline-rate parameter 22 ' 70,71 . 
The dashed line shows a linear fit to the data excluding SN 2004dt 22 . The colors of the 
symbols are the same as in Fig. 2 of the main text. 
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6 Supplementary Figure 3 
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The decline-rate parameter Am 15 (B) versus the late-time emission-line velocity shift. 
The three dashed lines schematically illustrate the expected viewing angle effect on the 
decline-rate parameter, which would produce a ~ 0.4 mag difference in the observed 
Am 15 (B) for an observer in the offset direction as opposed to the opposite direction 5 . 
The lines correspond to three hypothesized explosion configurations which are mutually 
different in M( 56 Ni) and therefore in the intrinsic luminosity. 
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